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ABSTRACT 

A  st  .’H;  is  made  of  the  influence  of  a  coplanar  mag¬ 
netic  field  on  the  stability  of  a  conducting  fluid  flowing 
between  parallel  planes.  After  derivation  of  the  general 
stability  equations  for  small  magnetic  Reynolds  number, 
numerical  results  are  obtained  for  the  case  where  initial  per¬ 
turbations  of  the  magnetic  field  vanish.  This  must  occur  if  the 
channel  walls  have  zero  or  infinite  conductivity.  Four  sets  of 
stability  diagrams  are  presented  so  that  each  stability  curve 
will  represent  the  effect  of  a  given  applied  magnetic  field 
as  only  one  of  the  four  quantities  in  the  Reynolds  number  is 
changed.  The  flow  is  always  stable  for  initial  disturbances  of 
the  field  produced  by  passage  of  a  pulsating  current  through  walls 
of  finite  conductivity. 


This  paper  was  presented  at  the  meeting  of  the  Fluid  Dynamics 
Division  of  the  American  Physical  Society  on  November  23-25,  1964 
in  Pasadena,  California, 
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I,  INTRODUCTION 

The  effei^t  oi'  ’J  uniform  cuplanai  magnet  in  field  on  the  stability 
of  parabolic  flow  of  a  conducting  fluid  between  parallel  wulls  has  been 
the  subject  cf  several  invest  iqat  i  oriu.  The  stability  equations  were 
first  fjiven  by  Michael  ^  ^  ^  who  showed  fjquire's  Thpcirem  is  .il-n  appilicabie 
to  two-dimensional  maqnetohydrodynamic  flows.  In  addition  to  the  Reynolds 
number  R,  tuo  other  nondimensiorial  {>ara»et»r-s  appear  in  the  equations; 
the  magnetic  Reynolds  riumber  R^,  and  the  interaction  parameter  N, 

Because  of  the  complicated  form  of  the  stability  equations, 

( 2) 

Stuart  simplified  the  equations  by  assuming  R^  is  small.  The  fluid 
dynamic  and  electrodynamic  equations  are  uncoupled  because  the  induced 
magnetic  field  is  of  second  order.  Stability  curves  were  obtained  for 
constant  values  of  the  parameter  q  =  to,  where  a  is  the  wave  number.  Un¬ 
fortunately,  the  drawbacks  of  this  form  cs  compared  to  the  use  of  N  alone 
were  not  recognized  by  Stuart  until  the  calculations  ware  well  under  way. 
The  curve  q  =  .08  formed  a  closed  loop  indicating  complete  stability  above 
a  certain  Reynolds  number,  Stuart  claimed  the  closure  of  the  curve  could 

not  be  verified  because  the  assumption  of  small  R  was  violated  at  the 

(3) 

higher  valuer  of  R.  In  a  discussion  of  Stuart's  paper.  Cowling  said 
the  existence  of  a  region  of  stability  for  large  values  of  R  "hardly  seems 
reasonable". 

Using  the  stability  equation  for  small  given  by  Stuart, 

( 4 ) 

Rossow'  '  solved  the  problem  again  and  obtained  stability  curves  for  ccn- 
stant  values  of  the  parameter  N,  The  curves  shift  to  higher  R  with  in¬ 
creased  N,  but  do  not  form  closed  loops.  This  suggested  the  closure  of 
Stuart's  curves  was  due  to  the  use  of  the  parameter  q  iriste  id  of  N,  This 
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u-as  shown  tc  be  the  case  by  Mains'  ,  who  obtained  quali*^ative  agreement 
between  the  rurues  of  Rossow  and  Stuart  by  crossplotting.  The  agreement 
is  only  goal itat iwe  because  Rossow  used  the  wroiu)  boundary  condition  along 
the  axis  of  the  channel.  The  boundary  condition  applicable  to  the  outer 
edije  of  a  flat  plate  boundary  layer  was  used  instead  of  tho  conditions 
for  an  ant  i  :.jymmetr  ic  disturbance  in  a  channel,  for  this  reason,  Rossow's 
results  are  not  considered  in  the  remainder  of  this  paper. 

Because  the  existing  stability  diagrams  are  of  limited  value, 
the  stability  equation  for  small  R^  has  been  resolved  by  an  exact  numeri¬ 
cal  method.  A  useful  set  of  stability  diagrams  plotted  for  several  non- 
dimensional  parameters  is  presented  in  this  paper.  In  addition,  this 
paper  attempts  to  clarify  several  points  which  have  led  to  some  confusion 
and  misunderstanding  about  the  effects  of  a  magnetic  field  on  stability. 

The  first  point  concerns  the  derivation  of  the  stability  equa¬ 
tions  for  small  R  ,  In  taking  the  limit  R  '♦0,  Stuart  neglected  a  term 
fli  in 

of  the  order  of  the  terms  retained,  and  arrived  at  an  equation  which  is 

specialized  for  the  case  where  initial  or  externally  produced  perturbations 

of  the  magnetic  field  vanish.  This  condition  is  always  satisfied  by  channel 

walls  that  are  perfect  conductors  or  insulators.  Because  B  is  small,  the 

m 

Field  influences  the  flow  to  first  order,  but  the  flow  influences  the  field 

only  to  second  order  so  that  first  order  perturbations  of  the  field  can 

only  be  generated  by  passage  of  an  oscillatory  current  through  channel 

walls  of  Finite  conductivity.  In  an  attempt  to  be  more  rigorous,  Tatsumi^^^ 
introduced  a  new  variable  but  cbtained  Scuart's  specialized  equation.  In 
effect,  he  also  neglected  first  order  perturbations  of  the  magnetic  field. 

In  this  paper,  the  general  stability  equation  for  small  R^  is  presented, 
and  solutions  are  given  for  channel  walls  with  zero,  infinite  and  finite 


conduct i V j  ty 


The  second  point  ccinrems  the  chctice  of  parameter  used  to 

describe  the  stability  curves,  and  the  correct  interpretation  uf  the 

curves  for  the  particulir  parameter  rhuseri.  In  suability  problems  where 

(7) 

the  Reynolds  number  is  not  the  only  nondimensional  parameter,  Hairs 
has  shown  the  charuje  in  H  should  be  ii.terpreted  as  a  variation  of  only 
those  quantities  not  common  to  both  R  arid  the  other  norid  imens  ional 
parameters.  For  the  parameter  q  chosen  by  Stuart,  this  would  mean  varia- 

A,  , 

tion  of  the  coefficient  of  viscosity  p  to  channe  R.  Since  R  is  inde- 

m 

pendent  of  u,  the  assumption  of  small  R^  is  not  v ini  ited  by  variation  of 
R  in  this  manner,  and  the  closure  of  Stuart’s  curves  is  therefore  correct 
if  R  =  R(u).  If,  as  Stuart  assumed,  R  -  R(u^),  the  closure  of  the  sta¬ 
bility  curves  cannot  be  decided  with  the  present  theory  because  the 
assumption  of  small  R^^  is  violated  for  laiye  values  of  R,  Since  varia¬ 
tion  of  to  chance  R  would  also  change  q,  some  other  quantity  in  q  must 
0 

vary  along  each  of  Stuart's  stability  curves  in  order  to  keep  q  fixed. 
Clearly,  some  new  parameter  independent  of  u^  must  be  chosen  to  describe 
the  stability  curves  if  R  =  R(u^).  In  this  paper,  parameters  are  intro¬ 
duced  which  permit  variation  of  each  quantity  in  R  individually. 
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II.  STABILITY  EQUATIONS 

The  magnetohydradynarnic  equations  which  govern  the  flow  of  an 
incompressible,  viscous  fluid  are  (in  dimensional  form): 

flow  continuity:  V  ,  V  =  0  (1) 

DU  ^ 

Navier-Stokes:  p  ■“  +VpijXB  +  |2V\/  (2 

Dt 

Ohm's  Law:  J  =  'SlZ  ~  ( 


Faraday's  Law: 

V  X  E  =  _ 

dt 

magnetic  continuity: 

7.0  =  0 

(b 

Ampere's  Law: 

-•  -* 

J  =  V  X  B 

(c 

where  p  and  a  are  the  viscosity  coefficient  and  electrical  conductivity 
respectively*  ^  is  the  fluid  pressure,  ^  is  the  density,  T  is  the  time 
and  C  is  the  electric  field*  The  components  of  the  magnetic  field  and 
velocity  are,  respectively: 

-  -  ^  *  W®2 

V  =  iu  ♦  jv  (P 

Elimination  of  C  from  Eqs*  (3)  and  (4)  with  the  aid  of  (1),  ( 
and  (6)  yields 

^  3  1  7^0  +  B»VV-\/*7B  (9 

dt  d 


A  second  equation  is  obtained  by  elimination  of  J  from  Eqs*  (2)  and  (6) 


P  iil  +  Vp  =  (B  .  V)B  -  '^B  ♦  4^  V  (10) 

Dt 


It  is  convenient  to  nondimencionalize  by  definini)  the  ntu-  v..riables 


where  L  is  the  channel  half-width,  and  u  is  the  velocity  at  the  axis  of 

o 

the  channel.  The  tilda  ('-)  indicates  a  dimensional  quantity, 

Cqs,  (9)  and  (10)  are  linearized  by  assuminq  small  perturbations 
in  the  dtapendeht  variable  of  the  form 


u  =  U(Y)  4  « 

dy 

dK 

p  *  Ptxty)  ♦  «a'(K,y,t) 
B,  =  1  ♦  c 

1  dy 

iliajxai 

2  dx 


where  e  is  a  small  quantity.  These  relations  automatically  satisfy  the 
continuity  Eqs,  (1)  and  (5),  Assuming  solutions  of  the  .eparable  form 


d(y)exp  ia(x  -  ct) 


(13) 


♦  z 

=  ♦(v)exp  ia(x  -  ct)  (14) 

The  linearized  form  of  Eqs,  (9)  and  (10)  reduce  to 

-  a^>ir1  =  0-(U-c)>|r  (15) 

aR  L  yy  J 

m  " 

_ D0  =  laNr(U-c)\|f  -  0] _ (16) 

where  0  is  the  operator  defined  by 

00  =  (U-c)(0  -a?0)  -  U  0  ♦  r0  -  2a^0  ♦  01^0!  (iV) 

yy  yy  “R  -  yyyy  yy  -1 

2 

The  steady-state  velocity  profile  is  U  >  1  -  y  •  These  two  coupled  equa¬ 
tions  are  the  stability  equations  for  plane  I*I1D  flow  with  a  uniform  co- 
planar  magnetic  field.  In  the  remainder  of  this  paper  the  operator  D  is 
used  so  that  the  le*'t  hand  side  of  Cq.  (16)  is  expressed  simply  as  00. 

The  boundary  conditions  along  the  channel  walls  require  the 
velocity  components  to  vanish 

- - 445)^ 

If  the  walls  are  perfect  conductors, 

♦(±  1)  =  0  (19) 

and  if  the  walls  are  perfect  insulators, 


t(±  •) 


0 


(20) 
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Tor  walls  oF  finite  conductivity,  *.he  value  of  1)  taxes  on  some  unknown 
value  found  by  matching  the  solution  in  the  channel  to  the  solution  in  the 
wall. 

Because  of  the  complicated  form  of  Cqs.  (15)  and  (16).  we  re¬ 
strict  our  attention  to  the  special  case  of  small  R  .  Assuming 


♦ 


♦ 


R 

m 


(21) 


<t  ♦  R  i  ,♦  ••• 
^3 - m  1 


Cqs.  (15)  and  (16)  become 


♦  -  = 

oyy  o 


0 


(22) 

(23) 


ThJi  is  the  general  form  of  the  stability  equations  for  small  R^. 

(2) 

Stuart  '  obtained  a  single  equation  which  corresponds  to  Eq. 
(23)  without  the  term.  It  la  therefore  specialized  for  the  case  where 
the  perturbations  in  the  megnetic  field  of  the  order  of  vanish.  In¬ 
stead  of  using  Cq.  (15),  Stuart  performed  an  order  of  magnitude  analysis 
on  the  equation 

0 

•  #  d  -  (U-c)t  (24) 

m  “o 

where  i  and  d  s  u^0  are  dimensional  quantities.  By  assuming  the 

terms  t  arid  0  instead  of  t  and  0  ware  of  the  same  order  he  neglected  the 

>1^  term  on  the  right  hand  side  of  Eq.  (23)  in  comparison  with  the  other 

terms.  Eqs.  (22)  and  (23),  obtained  by  assuming  ^  and  0  are  of  the  seme 

o  o 

order,  reduce  to  Stuart's  stability  equation  when  =  0. 
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In  an  attempt  to  be  more  rigorous  in  the  derivation  of  the  eta* 
bility  equations,  Tatsumi^^^  introduced  a  new  variable 

get  -  (25) 

I  yy 


and  used  tq.  (15)  in  the  form 


0f|  ♦  iaNti  =  ia(U-c)R|^Dt  (27) 

For  small  R  ,  Tatsumi  neglected  the  right  hand  side  of  Fq.  (27)  and  ob- 

R) 

tained  an  equation  of  the  form  of  Cq.  (29)  below.  He  proceeded  to  outlio 
the  solution  of  the  differential  equation  for  the  boundary  conditions 
q  =  q'  s  0  at  both  boundaries.  It  is  easy  to  show,  in  the  limit  as  R^  ' 
that  Tatsumi 'e  equation  is  satisfied  by  q  «  0.  By  using  Eq.  (21)  and  a 
similnt  power  series  sxpension  for  q,  Cqs.  (26)  and  (27)  are,  to  firct 
order, 

q  s  0  (2B) 

'o 

Oq^  ♦  ° 

The  corresponding  second  order  equations  ars 

iq„ 

0  .  —a  ♦  (u  -  c)+  (30: 

O  Qs  O 

Oq^  ♦  iaNq^  =  io(U-c)Di^  (31. 

It  is  clear  that  Tatsumi's  equation  is  not  equivalent  to  Stuart's  stabil 


Much  of  tho  con« 


equatior  because  the  variable  n  vanishes  for  ^cro  R  • 

*  IB 

fusion  results  because  Eq.  (29)  and  Eq.  (31)  with  =  0  have  the  same 
form  as  Eq.  (23)  with  s  0.  These  equations  are  by  no  means  equivalent, 
Tatsumi's  equation  is  meaningful  only  for  small  values  of  R^  greater  than 
zero,  while  Stuart's  equation  is  valid  for  R^  =  0  as  well.  Both  stability 
equations  assume  =  0,  which  is  a  special  case. 
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III.  CHANNEL  WITH  MALLS  OF  FINITE  CONDUCTIVITY 

If  the  walls  ha\/e  finite  conductivity  or^,  initial  perturbations 
of  the  magnetic  field  can  be  induced  by  passage  of  a  pulsating  current 
through  the  walls.  In  the  walls,  which  extend  to  y  =  1  the  magnetic 

field  must  satisfy  Eg.  (9)  with  V  =  0.  Using  Eg.  (12)  for  the  field  com¬ 
ponents  and  B^,  and  Eg.  (14)  for  the  form  of  the  disturbance,  Eg.  (9) 
reduces  to  the  form 

♦  ♦  ♦  =  0  (32) 

oyy  o  ' 

where 

♦  iaKc  (33) 

This  eguation  applies  only  when  R  is  small.  In  addition,  the  nondimension 

Rl 

parameter  K  s  9^  f  is  assumed  to  ba  of  order  one. 

If  the  disturbances  are  symmetric  with  respect  to  the  channel 
axis,  the  solutions  of  Cqs.  (2t)  snd  (3?)  are,  respectively. 


O^y^l  t«A  Cosh  ay  (34) 

o 

- y-i-4 _ AXoa  fty  _  (  j5) 


The  boundary  condition  requiring  the  matching  of  and  at  y  =  i;  1 
leads  to  the  relation 

a  tanh  a  «  -  P  tan  P  (36) 

Since  a  is  real  and  positive,  P  must  also  ba  real.  Eq.  (33)  shows  c  must 
be  a  negative  imaginary  number  since  K  is  positive.  The  disturbances  are 
always  stable. 
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l\l,  CHANNEL  lUITH  U/ALLS  OF  ZERO  OR  IMFINITE  CONDUCTIVITY 

If  the  channel  walls  are  perfect  conductars  =  °°)  or  perfect 

insulators  (a  =  O),  initial  perturbations  in  the  magnetic  field  must 
w 

v/anish.  vk  can  also  vanish  for  walls  of  finite  conductivity  if  the  current 
o 

J  is  zero.  The  eigenvalue  problem  reduces  to  the  solution  of  Eq.  (23) 

with  'I'  =  0.  The  boundary  conditions  are  given  by  Eq.  (IB).  As  is  cus- 

o 

tomary  in  the  nonmagnetic  case,  we  restrict  our  analysis  to  antisymmetr ical 
disturbances  so  that  0  is  an  even  function  of  y.  Integration  is  only  neces¬ 
sary  over  half  the  channel  if  the  boundary  condition 

0=0=0  (37) 

y  yyy 

is  satisfied  at  the  channel  axis. 

The  problem,  as  it  is  presently  posed,  was  first  solved  by 

Stuart^^^  using  the  Tollmien-Schlichting  theory.  In  this  paper,  solutions 

( 8 ) 

are  obtained  by  an  exact  numerical  method  developed  by  Hains  and  Price^ 
for  Poiseuille  flow  between  flexible  walls.  Details  of  the  numerical 
method  can  be  found  in  Reference  8,  but  a  brief  outline  will  be  given  here. 

The  numerical  solution  is  begun  by  dividing  the  channel  half-width 
into  n  parts  of  equal  length.  After  introduction  of  a  new  dependent  vari¬ 
able  developed  to  reduce  truncation  errors,  the  differential  equation  and 
boundary  conditions  are  written  in  finite  difference  form  at  each  of  the 
divisions.  The  resultant  system  of  linear  algebraic  equations  possesses 
a  nontrivial  solution  if,  and  only  if,  the  determinant  of  the  coefficients 
vanishes.  This  condition  determines  the  eigenvalue  c.  for  each  combina¬ 
tion  of  a,  R,  N  specified,  three  initial  guesses  for  c  are  made,  and  the 
value  of  the  determinant  for  each  guess  is  calculated,  A  complex  quadratic 
is  passed  through  these  three  values,  and  the  root  of  the  quadratic  nearest 
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the  last  guess  is  taken  as  the  new  approximation  to  the  elgenv/alue.  For 
n  =  50,  iteration  was  continued  until  the  value  of  the  determinant  had 
been  reduced  to  some  preselected  value  close  to  zero.  The  final  value  of 
c  is  the  eigenvalue.  This  numerical  method  was  programmed  for  an  IBM  7090 
computer  using  single  precision  with  8  digits.  The  values  of  c  given  in 
this  paper  are  accurate  to  at  least  four  decimal  places.  This  accuracy 
was  checked  by  repeating  the  calculation  of  some  points  with  a  double  pre¬ 
cision  pregram  carrying  16  digits, 

Meutral  stability  curves  for  the  wave  number  a  and  wave  propa¬ 
gation  velocity  are  shown  in  Figs.  1-6.  Three  sets  of  curves  are  pre¬ 
sented  for  the  three  nondimensional  parameters  l\I,  l\l  *  and  N*,  where 


N  -  NR  =  cr  L/p 

N'i'  H  n/R  =  a  p  By(p  u) 


I  (30) 

/ 


Since  N  is  independent  of  p,  these  parameters  have  been  selected  so  that 

N'  is  independent  of  P  and  u^,  and  N’*-'  is  independent  of  L.  Use  of  all 

three  parameters  to  describe  the  stability  curves  permits  the  variation 

of  R  to  be  interpreted  as  a  change  in  one  of  the  quantities  p,  u^,  L  or  p. 

In  Figures  1  and  2  where  N  has  been  chosen  as  the  interaction 

parameter,  a  given  value  of  the  applied  magnetic  field  will  lead  to  a  fixed 

value  of  N  if  the  viscosity  is  varied  to  change  R,  This  is  indicated  by 

R  =  R(p)  along  the  abscissa  of  the  figure.  As  the  figures  clearly  show, 

the  critical  Reynolds  number  increases  with  N  and  the  wave  propagation 

velocity  c  decreases  with  N, 
r 

In  a  channel  flow  experiment,  it  Ls  customary  to  increase  the 
Re/iiolds  number  by  increasing  the  flow  velocity  u  .  Stuart  also  interpreted 


Fig.  2  Neutral  stability  curves  For  the  wave  propagation  spaed 


Fig.  3  Neutral  stability  curves  for  the  wave  number  a  for  various 


0.26 
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Fig.  4  l\leutral  stability  curves  for  the  wave  propagation  speed 
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his  stability  curves  in  the  same  manner.  Instead  of  usiny  N,  or  to  as 

Stuart  used  to  describe  the  curves,  the  parameter  was  chosen  because 

of  its  independence  of  u^.  The  stability  curves  are  shown  in  Figs.  3  and 

4.  The  curves  are  similar  to  those  for  W,  but  approach  the  nonmagnetic 

curve  =  0  more  rapidly  with  increased  R(p,  Numerical  values  of  c 

1/3 

over  a  range  in  a  and  R  are  given  in  Table  I  for  some  representative 
values  of  N\  Positive  values  of  the  imaginary  part  of  c  indicate  an  un¬ 
stable  disturbance. 

The  last  set  of  curves  shown  in  Figs.  5  and  6  show  the  effect 
of  variation  of  the  channel  width  L  to  change  the  F?eynolds  number.  As  N* 
is  increased,  the  neutral  stability  curves  form  closed  loops  with  the 
region  of  instability  confined  to  the  interior  of  the  loop.  This  means 
the  flow  is  stable  for  small  channel  widths  and  for  large  channel  widths. 
When  N*  >  3,26  10  the  region  of  instability  disappears  and  the  flow 

is  stable  for  all  channel  widths. 

The  distribution  of  the  eigenfunction  is  shown  in  Figure  7.  Only 
the  real  part  of  0  and  its  first  derivative  are  shown  because  the  imaginary 
parts  are  small  in  comparison.  The  magnetic  field  tends  to  reduce  the 
gradients  of  the  velocity  perturbations  in  the  region  near  the  wall.  Using 
the  eigenfunction  for  N  =  0,  the  streamline  distribution  shown  in  Figure  8 
was  calculated.  For  clarity,  the  bending  of  the  streamlines  was  greatly 
exaggerated  by  taking  €  =  1,  The  antisymmetric  form  of  the  disturbance  is 
clearly  evident  from  this  figure. 
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V.  CONCLUSIOMS 

A  study  has  been  made  of  the  stability  of  a  conducting  fluid  be¬ 
tween  two  parallel  planes  with  an  aligned  magnetic  field.  When  the  magnetic 
Reynolds  number  is  small,  three  different  methods  of  producing  a  disturbance 
are  possible.  The  initial  disturbance  can  be  in  the  fluid,  in  the  magnetic 
field,  or  in  both  simultaneously. 

When  the  channel  walls  have  finite  conductivity,  a  disturbance 
of  the  magnetic  field  can  be  introduced  by  passage  of  an  oscillatory  current 
through  the  walls.  This  in  turn  causes  the  fluid  to  oscillate,  but  this 
type  of  disturbance  is  always  stable. 

Neutral  stability  curves  have  been  presented  for  initial  distur¬ 
bances  produced  only  in  the  fluid.  The  effect  of  the  magnetic  field  is  to 
increase  stability,  but  the  particular  shape  of  the  stability  curves 
depends  on  how  the  Reynolds  number  is  varied.  Four  sets  of  stability  curves 
were  presented,  each  corresponding  to  a  variation  of  one  of  the  four  quan¬ 
tities  in  the  Reynolds  number.  The  stability  curves  form  closed  loops  when 
the  channel  width  is  varied  to  change  the  Reynolds  number.  The  loops  dis¬ 
appear  and  complete  stability  is  obtained  when  N*  ^  3.26  •  lO”^.  When  other 
quantities  are  varied  to  increase  the  Reynolds  number,  the  stability  curves 
are  shifted  to  higher  Reynolds  number  as  the  interaction  parameter  is  in¬ 


creased. 
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